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ABSTRACT 


Moored current meter observations were made along the 
central California continental slope from Pt. Piedras Blancas 
to the Farallon Islands. The study area covered an alongshore 
distance of approximately 290 km and examined the data from 
off Pt. Piedras Blancas (P4), Pt. Sur (P2 and P3), Monterey 
Bay (MB1 and MB2) and the Farallon Islands (D and E). Time 
and frequency domain analyses were performed on three time 
segments that included data from combinations of the above 
mooring locations based on common time periods and depths (350 
me- 500 m). Segment 1 revealed a mean poleward flow 
attributed to the California Undercurrent at all moorings. 
Segment 2 had a similar mean poleward flow but also exhibited 
an equatorward reversal at periods between 19.5 and 58.5 days 
at P2 and MB2 that was less apparent at P3 further offshore 
and appeared to be coastally trapped. Observed wavelengths 
were compared with simple wave theory indicating that the 
observed signal may have resulted from a combination of two 
theories i.e. a coastally trapped wave. The complex bottom 
topography between Pt. Sur and the Monterey Bay prohibited 
more rigorous comparison between theory and the observations. 
Satellite sea surface temperature imagery during Segment 3 
showed an anti-cyclonic meander with a 65 km radius outside of 


the Monterey Bay. Current and temperature records at P2 and 


nog 


MB2 indicated that the surface feature's position varied and 
influenced the currents at depths down to 500 m. Time domain 
empirical orthogonal functions were calculated for alongshore 
and across-shore components separately. The first two 
empirical modes accounted for between 81% and 86% of the 
alongshore variance in all segments and were attributed to the 
mean California Undercurrent influence (mode 1) and to the 
deviations from the mean state (mode 2). The first two 
across-shore modes explained between 66% and 81% of the 
variance and consistently demonstrated opposing current 
structures between the Pt. Sur and Monterey Bay moorings, 


likely due to the Monterey Bay meander influence. 
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I. INTRODUCTION 


A. EXPERIMENT BACKGROUND 

The paper’s objective is to examine and describe the 
alongshore variability of sub-thermocline currents (350 m - 
500 m) over the California continental slope. The 
investigation employs conventional frequency domain (spectral) 
analysis as well as time domain empirical orthogonal function 
(EOF) analysis to provide a detailed picture of three 
different time segments. The analysis incorporated current 
meter data from several moorings between Point Sur and the 
Farallon Islands along the California coast (Figure 1). The 
study encompassed an alongshore distance of approximately 290 
km. All moorings were deployed on the middle slope between 
the 800 m and 2000 m isobaths. Mooring P4, the furthest south 
of all moorings was located off Point Piedras Blancas, 
approximately 100 km to the south of Point Sur. Two moorings, 
P2 and P3 were positioned off Point Sur. Mooring P3, provided 
data from the lower continental slope (1800 m) while P2, 
approximately 26 km east of P3, was positioned over the upper 
Slope (800 m). MB2, located 41 km to the north of P2, was 
positioned outside of the Monterey Bay (1979 m) while MB1, 
approximately 31 km to the east of MB2 was inside the Monterey 


Bay over the axis of the Monterey Canyon. Positioned 
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Figure 1. Current Meter Mooring Locations: A map of the 


central California coastal region. Moorings are 
shown as solid circles with alphanumeric labels. 


approximately 111 km to the north of MB2, mooring D marked the 
southern boundary of the Farallon Islands data set. Finally, 
mooring E was located 32 km north of D in the vicinity of the 
Farallon Islands outside San Francisco Bay. (Figure 1) 

The analysis of alongshore propagation required common 
time segments from each location. Various organizations 
collected data for different purposes, limiting the number of 
common time segments available. Instruments at depths between 
350 m and 500 m provided the most continuous alongshore data 
records so the analysis focused on these depths. The 
experiment used three segments of 110, 199 and 177 days that 
spanned an overall time period from December 1989 to November 


1991 (Figure 2). 
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Figure 2. Current Meter Time Line: Common time periods are 
labeled as Segments 1, 2 and 3 and are shown by 
different cross-hatching styles. 


The remainder of this chapter provides background 
information on the California Current System and describes 
various deviations from the mean flow observed there. A 
discussion on coastal trapped waves is presented here to 
include the derivation of the dispersion relationship for 
internal Kelvin waves. The second and third chapters review 
Pata collection and analysis techniques respectively. The 


fourth chapter presents the results and is followed by a 


discussion chapter that applies theory to the observations. 
The final chapter gives the conclusions with recommendations 


for future research. 


B. THE CALIFORNIA CURRENT SYSTEM 
1. The Mean State 

The data were collected from within the California 
Current System (CCS), the eastern boundary current of the 
North Pacific gyre. This sub-tropical gyre consists of an 
anti-cyclonic surface flow bounded by the westward flowing 
North Equatorial Current (20% N) to the south, the poleward 
flowing Kuroshio current to the west and the West Wind Drift 
to. Ehe montes When the West Wind Drift strikes the North 
American continent, it splits to the north to feed the 
poleward flowing Alaskan current and to the south to form the 
California Current. (Wooster and Reid, 1963) 

Large scale atmospheric forcing creates this anti- 
cyclonic circulation. During the Coastal Ocean Dynamics 
Experiments (CODE 1 in 1981 and CODE 2 in 1982), Halliwell and 
Allen (1987) attribute persistent, summertime equatorward wind 
stress to the interaction of the North Pacific subtropical 
high with the persistent low over the southwestern United 
States. The pressure differences between these two stationary 
systems are largest during the summer and hence equatorward 
wind stress is also largest in summer (Huyer, 1983).  Huyer 


(1983) describes the  southwestward migration of the 


subtropical high to about 28% N in February. Halliwell and 
Allen (1987) attribute winter time poleward wind stress to 
northward propagating cyclones and anti-cyclones. In general, 
atmospheric forcing drives an equatorward surface flow during 
summer which is often described as slow, broad and shallow 
(Wooster and Reid 1963). When these equatorward winds relax, 
the current reverses within 150 km of the coast (Winant et 
al., 1987; Chelton et al., 1988). This winter reversal is 
known as the California Countercurrent or the Davidson Current 
north of Point Conception (Chelton, 1984). 

Since this paper focused on the 350 m to 500 m depth 
interval, the poleward flowing California Undercurrent (CUC) 
influenced the data set more than any other feature of the 
CCS. Hickey (1979) observed the CUC at a depth of 300 m off 
southern Washington in 1972.  Chelton (1984) also observed 
this flow off central California and found an annual variation 
at Point Sur, in contrast with the semi-annual variation in 
the northern CCS such as that observed by Hickey (1979). 
Wickham et al. (1987) observed a "jet-like core" with velocity 
greater than 15 cm/s confined to within 30 km of the coast 
near Point Sur. Like Chelton (1984), they also found that the 
CUC varied annually with maximum poleward currents during May 
amar June of 1979 and' 1980. Chelton et al. (1988) found 
poleward flowing currents from the surface down to 400 m as 
far as 100 km beyond the central California shelf break during 


the summer of 1981 and 1984. They attributed this unusual 


occurrence to an extended period of anomalously weak 
summertime equatorward winds. Such poleward flow had been 
previously observed only during winter as the Davidson 
Current. 

Hickey (1989) studied the shelf and slope undercurrent 
off Washington and the Southern California Bight, nearly 
encompassing the entire CCS. Of particular relevance to this 
paper, Hickey discussed interruptions or anomalous cases where 
the generally poleward flowing undercurrent became 
equatorward. For example, in winter off Washington, 
alongshore wind stress and pressure gradient were opposite in 
direction and reversed from the summer case creating a 
wintertime equatorward undercurrent off the Washington slope. 
Additionally, Huyer et al., (1984) related equatorward flow 
over the slope to the presence of an eddy. Equatorward 
deviations likewise occurred in this data set and will be 
discussed on a case by case basis. 

2. The Baroclinic Equatorward Jet 

Contrary to Wooster and Reid’s (1963) California 
surface current description, data from the Coastal Transition 
Zone (CTZ) experiment revealed a strong baroclinic jet (core 
speed > 50 cm S^!) between 37? N and 38? N (Ramp et al., 1991; 
Huyer et al., 1991; Kosro et al., 1991). Huyer et al. (1991) 
superimposed acoustic doppler current profiler (ADCP) vectors 


onto satellite imagery and clearly showed a highly energetic 


meandering surface flow. This jet displaced isopycnals to 
depths greater than 500 m (Ramp et al., 1991). Daggett (1989) 
combined CTD and ADCP data to describe an equatorward 
meandering jet off Point Arena, California during CTZ. 


1 was observed in 


Equatorward velocity in excess of 25 cm sec” 
a shallow (50 m) core with weaker equatorward flow extending 
BHownweto 500 m. Additionally, in May 1988 Tisch (1990) 
observed an equatorward jet with velocity exceeding 25 cm 
sec"? off Point Sur down to 200 m. This flow was attributed 
Bo strong (-2.0 dyne cm ^) wind stress present during the week 
prior to hydrographic observations. In November 1988, 
equatorward geostrophic velocity was observed from the surface 
to depths exceeding 600 m (Tisch 1990). This paper provides 
additional evidence that an onshore meandering of this 
equatorward jet can interrupt the otherwise persistent 
poleward flow at depths down to 500 m (Chapter IV). 
3. Low Frequency Current Fluctuations 

In addition to the baroclinic jet interaction, several 
researchers have observed more periodic current fluctuations 
at periods from 5 to 40 days. Kundu, Allen and Smith (1975) 
compared dynamic and empirical modes of the velocity field 
near the Oregon coast. Although their observations were 
limited to the shelf (100m), they found high mutual coherence 
between the wind stress and barotropic u and v modes at 


periods between 7 and 16 days. Battisti and Hickey (1984) 


found significant coherence between two synthetic sub-surface 
pressure (SSP) levels over the entire low frequency spectrum 
band for coastally trapped waves (CTW) with periods greater 
than 5 days. They calculated the offshore-vertical structure 
and dispersion relationship of a free CTW and compared model 
results with observations. Again, these results were largely 
confined to the surface flow. 

Romea and Allen (1983) used an f-plane model to study 
vertically propagating, Coastal trapped, internal Kelvin waves 
(IKW). They paid particular attention to low frequency 
behavior with periods between 20 and 40 days, using three 
different types of wind forcing: an impulse wind at a point, 
a steady wind at a point, and a travelling wind. Their model 
showed that poleward travelling winds forced vertically 
propagating IKW while equatorward travelling winds forced 
waves that became surface trapped. They addressed the 
circumstances under which waveS may propagate vertically 
without dissipation and their results will be compared to the 
observations in this paper in greater detail. Vertical IKW 
propagation must be considered here as it provides a viable 
means by which a remotely forced wave could propagate between 
Point Sur and the Monterey Bay. Denbo and Allen (1987), 
Chapman (1987) and Lopez and Clarke (1989) found that remote 
atmospheric forcing to the south was an important source of 


Current variability on the western U.S. continental shelf. 


C. COASTAL TRAPPED WAVE THEORY 

A coherent, low frequency fluctuation (T = 19-29 days) 
persisted during the second time segment of this paper's data 
set. The dispersion relationships for simple internal Kelvin 
waves and barotropic continental shelf waves are presented as 
possible points of comparison for these observations. 

By using the inviscid momentum equations with the 
continuity equation, a series of wave-like solutions whose 


amplitudes decay exponentially with distance from a straight 


coastline can be produced. These solutions are generally 
referred to as "coastal trapped waves" (CTW) (Huyer, 1990). 
Huyer (1990) described two extreme types of CTW. On one 


extreme is the continental shelf wave (CSW) which assumes a 
homogeneous fluid and relies solely on a sloping bottom for 
trapping. This wave is also referred to as a topographic 


Rossby wave and has the following dispersion relationship 


(Gall 1982): 
TM 2£kA 
GC EQ AS) 
where o is the wave frequency in radians Secus p^ qs Ehe 


Coriolis parameter, k and 1 are the across-shore and 


alongshore wave numbers respectively and A is the bottom 


Slope parameter = - a where x is the distance offshore and 


ES 
H 


H is the bottom depth. The derivation of this relationship 


assumes a depth which increases exponentially with distance 
from the coast. 

On the other extreme, Kelvin wave theory assumes a flat 
bottom with a vertical boundary to trap waves of sufficiently 
low frequencies 1.e., W<<«< f. Internal Kelvin waves (IKW) 
represent the higher order modes of the Kelvin wave solution 
and depend only on the fluid’s density stratification (Huyer, 
TIGO) The IKW dispersion relationship is derived below 
following Gill (19829: 

Applying a vertical boundary at x = 0, all u (across- 
shore) components equal zero (i.e., current may not flow into 


the wall) and the following governing equations result: 


ZF = _i Op (1) 
p Ox 

dv __1 8p ; 

Ot p oy dd 
9w . 190p  p'g (3) 
dt p oz p. 

dv Ow 

AA ZTE 4 

am * az 0 (4) 


10 


= 0 (5) 





Equations (1) through (3) represent u, v, and w momentum 
respectively while (4) and (5) represent continuity and the 
equation GE state respectively. They assume an 
incompressible, inviscid linear flow where x, y, z, u, v, and 
w refer to the alongshore, across-shore and vertical 
directions and velocity components respectively, p is 
pressure, p, is the mean density under the Boussinesg 


/ 
approximation, PF is the reduced gravity (i.e., assuming 
O 


water column is stratified), f is the Coriolis parameter, and 
t is time. 


A pressure response form similar to that of Battisti and 


Hickey (1984) 


p= P(z) e`™* eiHy-ot] (6) 
is assumed where w equals wave frequency in radians sec !. 


Solving for v from equation (2) yields: 


v = BIST aek all-a (7) 


WP, 





Combining equations (3) and (5) produces a second order, 


linear differential equation for w: 


D 


32 s NONE -iW Or a n (8) 
Qe p m dz 





| d 
Where N^ is the Brunt-Vaisala frequency equal to ies << 
Z 


O 


radians sec"*. Solving for w: 


io dp 


w= ——> 
Po (N?-%?) az 


(9) 


The v and w velocities are then eliminated by substituting 
(7) and (9) into (4) producing a second order differential 
equation in P: 


d¢p 2 1*(N?-o?) 
dz- W? 





p sz (10) 


The rigid lid and flat bottom boundary conditions state 
that w= 0 at z = 0 and -H (where H is the bottom depth). The 
solution to equation (10) takes a sinusoidal form and applying 
the above boundary conditions yields: 

2 2:5 2 2 
Q^ 
where n is the internal mode number. Finally, solving for o^: 
PRG 
oe eS eee 
P na (12) 


(Ry er 


This relationship and that for a CSW are used to examine the 
coherent low frequency current fluctuations observed during 


Segment 2 of the data set (Chapter V). 
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II. DATA COLLECTION 


A. CURRENT METER DATA 

Current data were collected at seven different sites from 
Pt. Piedras Blancas to the Farallon Islands. Three separate 
time segments were chosen based on common periods when the 


data overlapped. Tables 1 through 3 summarize these segments. 


TABLE 1. CURRENT METER MOORING LOCATIONS AND DATA 
SOURCES: SEGMENT 1 (17 DEC 89 - 06 APR 0) a 


Bottom | Meter Current | Rotation 
depth depth Meter Axis 
E q PE CT 
EM (2r 44.910 E NE NPs 
fez | 00 _| 350 _| 36 - zoom |122- to. | ems | see | ws | 
A E les Ve Tel 
es | 1800 | 3s0_| 36" zoom | raz" av.ow | rows | see | mes 
es | 1800 [so |e- zo.o | 122 - 27.6% | eme | 30 | ws 
m [eso [soo [ae - aera | zr - onam [mors | sez | meari 
















Data 
Source 


















Since current meter depths from 350 m to 500 m provided the 
most complete and continuous data set, the analysis focused on 
these depths. Although other current meters were deployed, 
only the instruments at the common depths of interest are 
shown. Most of the data were collected with the Aanderaa RCM 


8 current meter with exceptions noted in Tables 1-3. 


T3 


Table 2. CURRENT METER MOORING LOCATIONS AND DATA 
SOURCES: SEGMENT 2 (30 AUG 90 - 17 MAR 91). 
Meter axis 


pany a 
Type (FT) 


mr er Ime faeces [or esa fame [we [ws | 


1800 POSTAN 27.6'W RCM 8 


m fes e [cene zoe E a P 
1979 rins 











Data 
Source 


Mooring Current Rotation 
























Table 3. CURRENT METER MOORING LOCATIONS AND DATA 
SOURCES: . SEGMENT 5 3^ ZA MA A RS NOV Ee vt 


— — 


Mooring Latitude Longitude Rotation 

axis 

(T) 
pean | eoz | 400 [se - 20.17 | 122% - 10.300 | romeo | saa | ws 
mer (eso Pare | 36 -4u.30u | 1220 - o2.5ów | remo | sez | mar | 


me | ame [em |se- comun | 122 22.5 memo | anat 
La oo sr -nam | 123 - 16.00 fae cms m 


1975 400 308-59" N 123° - 17.4'W. |M RCM 8 





* [ndicates linear interpolation between 350 m and 500 m 


1. RCM 8 Velocity Measurements 
A large vane assembly (48.5 cm by 50 cm) kept the 
instrument aligned into the current. A magnetic compass 
inside the recording unit recorded the direction of 
orientation via a potentiometer ring fastened to the lid of 


the compass. Current direction was resolved to 0.35% with an 
accuracy of + 5° for current speeds from 5 to 100 cm/sec ands 7.5° 


for speeds from 2.5 to 5 and 100 to 200 cm/sec (Aanderaa, 


1987). Combined magnetic deviation and variation of each unit 
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was determined on a surveyed test bench at the Naval 
Postgraduate School prior to each deployment. 
A shrouded paddlewheel with magnetic coupling through the 
top end plate measured current speeds in a range from 2 to 250 
cm/sec with accuracies of + 1 cm/sec or + 2% of actual speed, 
whichever was greater. The instrument performed a vector 
averaging by sampling the number of rotor revolutions and 
direction every 36 seconds for a 30 minute sampling interval. 
The meter maintained and intermediately stored a running sum 
of North-South and East-West speed components until the data 
sampling interval elapsed. At this time, the resultant 
average speed and its angle were calculated and recorded 
(Aanderaa, 1987). 
2. Temperature Measurements 
The RCM 8 used a thermistor with a variable resistance 
(5700 ohm at 0% C and 2001 ohm at 25° C) for temperature 
measurements. The thermistor is connected in series with low 
range (-2.469 C to 21.48? C), high range (6.089 C to 31.04? C) 
and wide range (-0.34° C to 32.17° C) resistors. Temperature 
values were then determined by various polynomials’ for 
corresponding temperature ranges based on the measured 
resistances. Temperature readings were resolved to 0.1% of 


the range and accuracy was +0.06° C. 


T5 


3. RCM 5/RCM 8 Differences 

Current meter moorings MB1, MB2 (Pillsbury, et al., 
1992) and D deployed the Aanderaa RCM 5 current meter. An 
older version of the RCM 8, the RCM 5 differed in several 
important aspects. First, instead of the solid state data 
Storage unit (DSU) used in the RCM 8, the RCM 5 relied ona 
1/4" magnetic tape and tape reader to store and retrieve data. 
Also, the RCM 5 used a Savonius type rotor instead of the 
Shrouded paddlewheel used on the RCM 8, The shrouded 
paddlewheel was designed to eliminate the unwanted effect of 
pumping the Savonius rotor through vertical heave in the 
mooring. Additionally, the RCM 8 vane area was redesigned 
from that of the RCM 5 to decrease the response time for 
current direction shifts. Finally, the RCM 5 does not use 
vector averaging. It recorded speed continuously but used a 
Snapshot direction at the end of the sampling interval, 


introducing possible errors in current direction measurements. 
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III. ANALYSIS METHODS 


A. NPS DATA 
1. Data Retrieval, Calibration and Editing 

The Aanderaa RCM 8 used a solid state data storage 
unit (DSU) to record current speed and direction, temperature, 
pressure and, when so equipped, conductivity. The DSU had a 
maximum storage capacity of 65,530 ten bit words, which 
corresponds to seven months of data when using a sampling 
interval of 30 minutes. All instruments used a 30 minute 
Sampling interval except those at the Monterey Bay moorings 
and mooring D, which sampled data every hour. The RCM 8 
determined sampling time with an internal quartz clock. 
Within all operating temperatures, the clock accuracy was 
sec day”?. 

Upon instrument retrieval, the DSUs were removed and 
downloaded to an IBM compatible personal computer (PC) via the 
DSU reader 2995. The reader converted bit codes to ASCII 
coded format (Aanderaa, 1987) and stored the converted data on 
the PC. Next, a calibration program applied conductivity, 
pressure, speed, time and direction calibrations to convert 
the raw data to engineering units of current speed (cm/s), 
direction (°T), water temperature (°C), salinity (PPT) and 


pressure (decibars). At this point, velocity records were 
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truncated to eliminate unusable data entries i.e., data from 
instrument failures and zero data records obtained before and 
after instrument deployment. The records were plotted and 
visually inspected for obvious outliers. A despiking program 
then checked the data to ensure all parameters were within 
certain tolerance ratios. 
2. Filtering 

A filtering program (Denbo et al., 1984) first 
employed a low pass (LP) Cosine-Lanczos filter (half power 
period at 2.9 hours) to remove high frequency noise from the 
data. Signals with periods greater than 3.6 hours passed 
unattenuated and signals with periods shorter than 2 hours 
were removed entirely. The program applied filtering weights 
to successive 25 point segments. A Lagrangian interpolation 
subroutine shifted the start time of the filtered series to an 
even hour. This interpolation scheme employed a third order 
polynomial whose coefficients were based on the desired 
fractional distance between two successive data points. Only 
data with 30 minute sampling intervals were LP filtered. This 
included Farallon (mooring E), Point Sur and Point Piedras 
Blancas data. 

A low low pass (LLP) filter (half power period at 46.6 
hours) used the output from the LP interpolated filter to form 
six hourly records. This filter, also a Cosine-Lanczos 


type, applied filtering weights to 121 point segments to 
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remove diurnal and shorter period energy as well as any 
inertial period energy. Monterey Bay data had an hourly 
sampling interval already and thus only used LLP filtering. 
The data sets now contained despiked 6 hourly records with 
tidal and inertial influences removed (Figure 3). 
3. Bridging 

Gaps were present in various time series due to 
instrument failure or during mooring turnarounds. Gaps of 
four days or less were bridged using a special program 
described below. All discontinuities greater than a few days 
were not filled. The bridge program first determined the 
width of the gap between the two input time series. With 
Lagrangian interpolation, the program then shifted the start 
time of the second series so that it started exactly 6 hours 
after the last time of the gap. A prediction error filter 
(Anderson, 1974) based on data from 25 days on either side of 
the gap determined the data points within the gap. The 
program accomplished this by extending the first time series 
3/4 of the way across the gap forward in time and the second 
time series 3/4 of the way across the gap backwards in time. 
It melded the overlapping data using the prediction filter, 
stored the data separately, and finally combined all data to 
mI the gap. Tests showed that this procedure had no 


detrimental effects on the autospectra. The requirement for 
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time series from mooring E 400 m with LP (top) 
and LLP filtering (bottom). 
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a continuous time series for spectral and EOF analysis made 
this procedure particularly useful. 
4. Coordinate Axes Rotation 

The interpretation of geophysical time series is often 
more straightforward in an alongshore vs across-shore rather 
than an east-north coordinate system. Thus, coordinate axes 
rotation facilitated across-shore vs. alongshore analysis. 
This procedure first required primary axis determination 
(Tables 1-3). A computer program determined the lengths of 
the semi-major axis, semi-minor axis and the angle from the 
positive x axis to the semi-major axis (Kundu, 1976). The 


program used demeaned u and v to determine the angle from the 


positive x axis, 8: 


N 
i > (covariance) 
0 = 1/2 tan“! = (13) 
N 
x DPI) 
L 1=1 
where covariance - u, * vi, u and v are the component speeds 


in the x and y directions respectively and N is the number of 
data points in the record. Based on these results, data were 
rotated locally with the exception of mooring P3. Because of 
the limited amount of data at P3 and to provide a consistent 


axis in the vicinity of Point Sur, P3 was rotated according to 
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the results from P2 and MB2. The rotation computed new 


component speeds based on 


E 
iI 


rot = U-COS() + v sin (@) 


Vrot 


v,cos ($) - u,sin(d) 
where u,g. and v,g. are the rotated component speed snm 


and y direction respectively, ug and vg are the mM 


O 


component speeds and $ is the degrees counterclockwise from 
true north (4 = 90 D from equation (13) above). The 


rotation process concluded the preliminary data processing 


procedures. 


B. OTHER DATA 

Portions of the data were collected by organizations other 
than the Naval Postgraduate School (NPS). Specifically, data 
from the Monterey Bay moorings were generously supplied by Dr. 
C.H. Pilskaln of the Monterey Bay Aquarium Research Institute 
(MBARI). The data were received as calibrated, hourly 
records and were subsequently filtered to form 6-hourly 
records (Pillsbury et al., 1992). The data set provided a 
critical piece linking the Farallon Islands and Point Sur 
current meter  moorings. Additionally, the Science 
Applications International Corporation (SAIC) worked in 
conjunction with NPS on a project that described the general 
circulation in the vicinity of the Farallon Islands. SAIC 


conducted calibration, maintenance, deployment and data 
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retrieval of all current meters at mooring D. The data at 
mooring D were provided by M. Noble (USGS), who carried out 
additional quality control, filtering and axis rotation. For 
continuity, all Farallon data (moorings D and E) were rotated 


eo 3289 T based on these results. 


C. TIME DOMAIN ANALYSIS 
1. Time Series Plots 
Vector velocities, across-shore and  alongshore 
velocity components and temperature time series were plotted 
for the common time segments. To facilitate alongshore 
comparisons, stack plots were formed presenting a picture of 
the moorings in a south to north direction; i.e., data from 
the southern most mooring from each segment is on top, the 
northern mooring on the bottom and the others in between. 
These plots proved critical in identifying the various events 
on which to focus additional analysis techniques. 
2. Empirical Orthogonal Function Analysis 
Time domain empirical orthogonal function (EOF) 
analysis (Kundu et al., 1975; Beardsley et al., 1985; Winant 
et al., 1987) provided a compact description of the alongshore 
variability of the current structure. The technique found the 
orthogonal eigenvectors and eigenvalues of a covariance matrix 
containing the desired variables. The number of time series 
equalled the number of eigenvalues, with each eigenvalue 


having a corresponding eigenvector or eigenfunction. The sum 


2:3 


of eigenvalues equalled the total variance among all 
components of the covariance matrix and thus each eigenvalue 
represented a percentage of the total variance for its 
corresponding mode. Modes were "ranked" according to the 
percent variance contained, i.e., the eigenfunction with the 
greatest variance was called mode 1, the next most mode 2, and 
so on. 

The technique allows one to decompose some parameter 
into a few primary modes representing the essential 
variability. Although meteorologists first employed this 
technique, its applicability to ocean phenomena became quickly 
apparent. Kundu et al. (1975) decomposed current velocity 
fields off the Oregon coast into both dynamical and empirical 
modes. They found that the empirical modes did not depend on 
the accuracy of simplifying dynamical assumptions but only on 
the statistics of the data. They also found that empirical 
modes were discretely orthogonal while dynamical modes were 
not. Beardsley et al., (1985) used both time and frequency 
domains to describe current and temperature variability in the 
Nantucket Shoals Flux Experiment (NFSE79). Finally, Winant et 
al., (1987) described moored wind, temperature and current 
observations made during the Coastal Ocean Dynamics 
Experiments (CODE) 1 & 2. In all cases, these researchers 
described the spatial structure of multiple variables ina 


powerful, compact manner. 
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3. EOF Theory 


Following Kundu et al., (1975) a brief outline of the 
BOF expansion technique is given. For this example, let 
E) denote demeaned, alongshore velocity at time tk (k = 
EFE K) and at mooring location x» {1 = 1,...,N). The 


matrix of time series at different mooring locations can be 
decomposed into a modal amplitude and a corresponding 


eigenfunction so that 


N 
v.(x,) = V E,(x)6, n=1,...,N 
n-1 


where , defines the empirical orthogonal eigenfunction and 


Acne modal amplitude for each record in the time series. 
Determining the least squares error between the actual 


velocity and the summation of decomposed velocities shows that 


N 
Yo R(x,,%;)6,(*;) 7» Ab (xj), m=1,...,N 
i-1 
where 
A K 
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represents the real and symmetric covariance matrix, A, 


corresponds to each eigenvalue and $4, obeys the orthogonality 


Condition: 


N 
» Oa (Xi) 6,(x,) E oT 
171 


where 6,, = 1 form = n and O form * n. The eigenvalues, A, 


represent the time average energy in each mode i.e., the 
variance and the sum of the eigenvalues both equal the total 
energy. 

A simple MATLAB function calculated the covariance 
matrix for the desired variables by first removing the time 
mean from each series. The function then found the 


eigenvalues A,, eigenvectors $4, and modal amplitudes, E, as 


described above. 


D. FREQUENCY DOMAIN ANALYSIS 

Spectral analysis employed the Fourier transform to 
provide energy estimates as a function of frequency. A 
modified version of the Matrix Laboratory (MATLAB) function, 
SPECTRUM computed spectral characteristics of across-shore and 
alongshore series. To allow investigation of the maximum 
frequency range without statistical degradation, the data 
segments were divided into five, equal length, overlapping 
pieces. A Hanning window, applied to each piece separately, 


damped the effects of the Gibbs phenomenon resulting from the 
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discontinuity created Dy the truncation of an infinite series: 


The coefficients of the Hanning window were computed by 


win) = 00851 = cosí2x EN eee N 





ea 


where w represented Hanning window weights and N the number of 
data points in each piece. 

Since the piece lengths were not powers of two, the 
modified program applied a mixed-base fast Fourier Transform 


(FFT) to each piece. The transform used by the FFT is: 


Ner 
A la 
n=0 


where X defined the transformed time series, x the Hanning 
windowed time series and N the number of records in each piece 
(Little and Shure, 1992). The program then calculated the 
power spectral density, cross spectral density, complex 
transfer function, coherence and phase functions and 
accumulated results from each piece. The pieces were averaged 
together to provide the final spectral estimates with 
confidence limits. Data were also plotted in variance 
conserving form, obtained by normalizing each spectral 
estimate by the frequency bandwidth and multiplying it by 2.3. 
The results from this and all analysis techniques are 


presented in Chapter IV. 
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IV. RESULTS 


Each time segment presented its own unique results. 
Specifically, currents appeared to vary seasonally between 
each segment. Fluctuations could have been influenced by 
surface seasonal transitions or perhaps deeper current 
variations. This seasonality coupled with distinct, mesoscale 
deviations from the mean flow made it sensible to describe 


each segment separately. 


A. SEGMENT 1 

The first time segment included data from 17 December 1989 
through 04 April 1990 (110 days) (Figure 2). Data were 
collected from current meters at P4 off Pt. Piedras Blancas, 
P2 and P3 Off Pt. Sur and MB1 in Monterey Bay. Instruments at 
350 m and 500 m provided data at all moorings except MB1 (500 
m only). P4 was located approximately 100 km to the south of 
Pt. Sur and 47 km offshore along the 1100 m isobath. Located 
in 650 m of water, MB1 was 35 km to the north of Pt. Sur, 14 
km offshore and very close to the axis of the Monterey Canyon. 
P2 was positioned 24 km offshore and 13 km from the shelf 
break (200 m isobath). Figure 4 schematically illustrates 


mooring locations and separation distances. 
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Figure 4. Segment 1 Mooring Separation Distances: 
Alongshore and offshore distances are shown with 
the approximate relative position of the Monterey 
Canyon. The canyon axis is indicated by the 
darker shading. 
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1. Time Domain Analysis 
Basic time domain current statistics for Segment 1 
were calculated for the time period (Tables 4a and 4b) and are 
described for each mooring below. At P4, basic statistics 
showed a mean  alongshore (ALS) component speed of 


E 


approximately 1.5 cm sec” at both depths with a maximum of 


T 


14.4 cm sec” and a minimum of -10.3 cm sec’! (Table 4a). 


The across-shore (ACS) component mean (-1 cm sec"*) was 
offshore at both depths but is likely within the error of the 
axis rotation. Mean currents were primarily in the alongshore 
direction. The maximum offshore current was -11.9 cm sec + 

while the maximum onshore current was 5.9 cm sec’? (Table 
"ono From 17 December until 06 January, a vertical shear 
existed at P4 between 350 m and 500 m with poleward flow at 
350 m and equatorward flow at 500 m (Figures 5 and 6). Flow 
became equatorward at both depths from 06 to 16 January when 
it turned poleward until late winter. Around 08 March, the 
current shifted to equatorward and remained so until the end 


of the segment. The currents looked coherent vertically from 


06 January to the end of the segment. 
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Table 4a. TIME DOMAIN STATISTICS FOR SEGMENT 1 
2... (ALONGSHORE COMPONENT) — — .,. — — 


Moor ing Meter Mean 
depth | (cm/s) 
(m) 
































Standard 
Deviation 
(cm/s) 


Mode 1 
Eigenfunction 
(Eigenvalue = 
172.0 (cm/s)2 
for 66.9% of 
total variance) 


Mode 2 
Eigenfunction 
(Eigenvalue z 
47.8 (cm/s)2 
for 18.6% of 
total variance) 












Table 4b. TIME DOMAIN STATISTICS FOR SEGMENT 1 
. (ACROSS-SHORE COMPONENT) — — COMPONENT) 
for 52.8X of for 20.4% of 


Moor ing Meter 
depth e 
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total variance) total variance) 
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Standard 
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(cm/s) 


Mode 1 
Eigenfunction 
(Eigenvalue = 
45.22 (cm/s)2 17.51 (cm/s)? 
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At P2, the mean alongshore current at 350 m for the 


entire segment was 6.8 cm sec! with a maximum of 28.4 cm sec ! 


and minimum of -5.9 cm secl. The 500 m instrument exhibited 


a similar mean (3.9 cm sec 1) with a maximum of 26 cm sec"? 
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1 


cm sec’! and minimum of -6.7 cm sec”! (Table 4a). The across- 


shore component mean was offshore (-0.7 cm sec"*) at 350 m 


1 and a 


with a maximum offshore velocity of -13.7 cm sec” 
maximum onshore velocity of 7.4 cm secl. At 500 m, the mean 
was onshore (0.5 cm sec"*) (Table 4b). The small across-shore 
means are again due to the coordinate rotation, which 
Statistically minimized the ACS component. 

P2 exhibited much stronger poleward flow than P4 from 
the beginning of the segment until 01 February. On this date, 
the poleward current abruptly diminished, then shifted to 
equatorward for 1 day at 350 m and 10 days at 500 m (Figures 
5 and 6). As suddenly as the current changed, it shifted back 
to poleward and reached the segment's maximum of 28 cm sec"? 
at 350 m on 16 February. This short relaxation was not 
observed at the other moorings. The current remained poleward 
until 05 March when the current at both instrument depths 
oscillated with a period of about 10 days for the final month 
of the time segment. The records from 350 and 500 m looked 
vertically coherent the entire time. 

The most consistent poleward flow occurred at P3. 
With a mean alongshore current of about 4.3 cm sec? at 350 
m and 3.4 cm sec”? at 500 m, the flow remained steadily 
poleward and offshore (across-shore component mean -3.4 and 
-2.6 for 350 m and 500 m respectively) from the beginning of 
the segment until 17 March (Tables 4a and 4b). On 9 


February, the 350 m current reached a maximum of 17 cm aec'h. 
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This occurred about one week prior to the maximum at P2. From 
10 February until 17 March, the flow at P3 rapidly shifted to 
equatorward at both depths. Additional data (not shown) show 
that this event exceeded 1000 m depth. During this event, 
current vectors rotated counterclockwise so that the flow was 
directed offshore until 01 April at which time the current 


shifted to onshore until the end of the segment (Figures 5 


I 


and 6). The maximum equatorward flow exceeded 24 cm sec * at 
350 m and 17 cm sec"? at 500 m. 
With a mean alongshore current of 1.2 cm sec? , the 


flow at MB1 500 m demonstrated significantly different 
characteristics throughout the time segment (Table 4a). The 


D 


poleward currents reached a maximum of 9 cm sec ^ on 03 March 


and a maximum equatorward flow of -7 cm sec! on 20 March. 


rur equatorward bursts appeared with a period 


Weak (= 5 cm sec. 
of about 12 days.  Mooring MB1 did not display the dominant 
poleward current observed at the other moorings (Figure 6). 
MB1's location near the axis of the Monterey Canyon seemed to 
isolate the mooring from the general flow pattern offshore 
during this time segment. 

The alongshore and across-shore velocity components 
were separately decomposed into empirical modes using the 
method outlined in Chapter III. Beginning with the alongshore 


component, the first two modes accounted for 85.5% of the 


total variance. The first mode alongshore eigenfunction 
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showed all moorings with the same sign (Figure 7, top left). 
This represented the strong poleward current seen at moorings 
P2, P3 and P4. The second mode alongshore eigenfunction 
indicated an opposing structure between P2 and P3 with P4 
again resembling the flow at P3 (Figure 7, bottom left). 
This perhaps represents two dynamical regimes (discussed 
later); an inshore or upper slope regime (including mooring 
P2) and an offshore or lower slope regime (including moorings 
P3 and P4). Mooring MB1 made only a small contribution to 
these modes, due to its small amplitude and lack of 
correlation with the other moorings. 

For the across-shore component, the first two modes 
accounted for 72.4% of the total variance (Table 4b). The 
first mode eigenfunction indicated opposing flow between Pt. 
Sur moorings (P2 and P3) and the Pt. Piedras Blancas mooring 
(Figure 7, top right). The second mode eigenfunction however, 
revealed an opposing structure between P2 and P3 with P4 
resembling the flow at P3 (Figure 7, bottom right). This 
again Supported the idea that P3 and P4, both located over the 
outer slope were contained within a somehow different 
dynamical regime than P2 and MB1. 

Alongshore modal amplitude and component velocity time 
series were correlated to determine which mooring dominated 
the different eigenmodes (Table 5). The correlation 
coefficients showed a very high correlation between P3 and 


mode 1 (.94/.95 at 350m / 500m). The P3 component velocity 


By 


time series and mode 1 amplitude were plotted and showed 


strong similarities (Figure 8). This indicated that the 


Table 5. SEGMENT 1 CORRELATION COEFFICIENTS FOR FIRST 
AND SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 


€ K- ee — ——À 


P4 350 m 
P4 500 m 
P2 350 m 
P2 500 m 
P3 350 m 
P3 500 m 


MB1 500 m 





current at P3 contributed the most to the mode 1 variance. 
Further, mode 2 showed the highest correlation at P2 (not 
shown) suggesting that P2 contributed most to the mode 2 
variance. 

The across-shore modal amplitude time series were also 
computed and correlated with each of the across-shore 
component velocity time series (Table 5). Consistent with 
the relative magnitudes presented in Figure 7, the first mode 
amplitudes were highly correlated with the P3 time series at 
both 350 m and 500 m (.98 and .96 respectively). The P3 350 
m across-shore component time series was plotted with the mode 
1 amplitude time series clearly showing the  current's 
influence on the mode 1 amplitudes (Figure 8). For mode 2, 


P4 showed the highest correlation (.71 and .62 for 350 m and 


38 


Segiment 1:49 $00 vs Mude 1 Amplitade (Alongshore Components) 


4n «e mre umi ami ee PR ee = oe a ere s UL 





Mode 1 


10 À 


fa | la Jue 


a E po l) ts : ? ANA | 
K A nal T ee! M dd " 
a x K i s 
11500 MI 
n EN 


\ 
S 
Ww A e. Ar 
fe » y 1 tt " r 14 ' , e 1 L 
mec pat TT un ira 
Segient 1 14 350 19 Mide Y Amplitude (Across shore Components) 
E nca DONE I lat M cr. 





20 [EET] veros 11T po. eee t 100109009? 
i Ll nm 


Figure 8. A Comparison of the Mode 1 Velocity Amplitudes 
with the Best Correlated Individual Time Serles 
for Segment 1: The top panel shows the 
alongshore components and the bottom panel shows 
the across-shore components. 


500 m, respectively) with flow at P2 opposing that at P4 and 


slightly less correlated. As expected, the flow at MB1 
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contributed very little to the overall structure in the first 
two modes with correlation coefficients close to zero 
(Table 5). 

To summarize, time domain analysis described a current 
structure characterized by primarily poleward alongshore flow 
with a small offshore component at P2 and P3 and a small 
onshore component at P4. A secondary influence caused 
opposing flow between P2 and P3 with P4 generally following 
the P3 pattern. An interesting result showed that although P2 
had a larger mean signature, P3 currents dominated the first 
empirical mode. The significance lies in the EOF analysis 
method where contributions from all moorings were considered 
to produce an overall spatial structure and will be discussed 
in greater detail later. Finally, results from MB1 suggested 
that this mooring was secluded from the offshore flow during 
this time segment. 

2. Frequency Domain Analysis 

Variance conserving spectra were plotted to determine 
the frequencies of maximum energy (Figures 9a through 9d). 
Since the currents were initially rotated in a manner that 
maximized alongshore variance, the alongshore component 
spectra generally showed the highest energy levels and are 
discussed first. Reported energy levels reflect the 


normalized values to be consistent with the plots. 
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Segment 1 Variance Conserving Auto-Spectra for 
spectra are shown for P4 
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Figure 9b. Segment 1 Variance Conserving Auto-Spectra for 
350m: Alongshore and across-shore components 
are shown for mooring P3 350 m. 
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Alongshore and across-shore component 
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Segment 1 Variance Conserving Auto-Spectra for 
spectra are shown for P4 


500 m: 
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Figure 9c. 
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Figure 9d. 
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Segment 1 Variance Conserving Auto-Spectra 


Across-shore Components 


Alongshore Components 


(500 m): Alongshore and across-shore component 
spectra are shown for P3 (top) and MB1 (bottom). 
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At P4, the 350 m current spectrum showed an energy 


2 


peak of 13.0 cm Secu e days (Figure 9a, top left). At 


500 m, the energy decreased and shifted towards shorter 


periods, peaking at 8.6 cm? S 


Sec with a period of 12 days 


igure 9c, top left). For P2, the maximum alongshore 
component energy level at 350 m was 28.3 cm? sec? and 


occurred at a period of 18 days . A secondary signal of 15.8 


2 22 


cm^ sec appeared at a period of 6 days (Figure 9a, bottom 


left). Similar peaks occurred at 500 m except the energy 


level at 6 days (24.4 cm? sec?) actually exceeded the 18 day 


2 2 


signal of 23.6 cm sec'^ (Figure 9c, bottom left). At P3 350 


m, the energy level decreased from that at P2 to a peak level 


2 2 


of 9.4 cm sec ^ at a period of 6 days (Figure 9b, left). The 


500 m spectrum here also exhibited a maximum energy level at 


2 2 


a 6 day period but decreased to 5.0 cm sec'^ (Figure 9d, top 


left). The MB1 alongshore spectrum (Figure 9d, bottom left) 


2 2 


displayed a maximum energy level of 6.4 cm^ sec ^ at a period 


Of 12 days. 
The across-shore component spectra exhibited energy 
peaks at periods similar to the alongshore components. The P4 


350 m spectrum peaked at 18 days with an energy level of 14.1 


2 =2 


cm“ sec (Figure 9a, top right). The energy maximum 


2 -2 


decreased to 4.4 cm“ sec at 500 m and occurred at 12 days 


üsggure 9c top right). At mooring P2, the 350 m spectrum 
showed two peaks again at 12 days and 6 days. The first had 


2 


an energy level of 8.7 cm sec"? while the second was slightly 
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z 2 


less energetic at 8.0 cm sec ^ (Figure 9a, bottom right). At 


500 m, the energy also peaked at 12 days but decreased to 6.1 


cm^ sec'^. A smaller 3.6 day peak occurred with an energy 


2 


level of 5.0 cm* sec". At P3, the 350 m spectrum peaked at 


9 days with an energy level of 9.2 cm* sec’? (Figure 9b, 
right). At 500 m, three distinct peaks at 9, 5, and 3.3 days 
occurred with energy levels close to 4.0 cm? sec"? at all 
periods (Figure 9d, "top rght s The spectrum at MB1 
displayed peaks at 4.5 and 3 days with energy levels around 
2.9 cm? sec'^ (Figure 9d, bottom right). Im generale 
across-shore component spectra at moorings P4 and P3 displayed 
energy levels comparable to those for the alongshore 
components while P2 and MB1 exhibited higher energy level in 
the alongshore component again suggesting two distinct 
dynamical regimes. 

In short, the alongshore component energy levels 
decreased with depth at moorings P4, P2 and P3 with the 
greatest energy at P2. The energy level at MB1 was comparable 
to the 500 m spectra at both P4 and P3. 

Cross-spectra between moorings showed relatively few 
coherent signals. For the alongshore components (Figures 10a 
and 10b) P2 and P3 were coherent at periods of around 5-6 
days. The signal at P3 led that at P2 by about 40? at 350 m 
and 169? at 500 m. The analysis also showed a coherent 4.5 - 


6 day signal between P4 and P2 at 500 m and a 6 day signal 


between P4 and P2 also at 500 m. P4 led P2 by 26%, 35? and 
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46° for signals of 6 days, 5.1 days and 4.5 days 
respectively. P4 led P3 by 13° for the 6 day signal there. 
The Significance of this 6 day fluctuation is not clear since 
it represented such a small fraction of the energy at P4. 
This may be a spurious coherence resulting from the strong 6 
day peak at mooring P2. No Significant coherence was observed 
at the dominant periods. 

The only across-shore coherent Signals were observed 
between P2 and P3 at both 350 m and 500 m. The 350 m data 
(not shown) exhibited coherence at 12 days, 5 days, and 3 
days. At the 12 day and 3 day signals, P2 led P3 by 5° and 
18° respectively while P3 led P2 by 29? at the 5 day signal. 
Coherence values just barely reached the 95$ significance 
level at this depth. Although the P2 energy reached a maximum 
at 12 days, the P3 energy at 12 days represented a much 
smaller fraction of the total variance; likewise for the 5 day 
and 3 day signals. Thus the importance of these coherent 
Signals is not clear. At 500 m, coherence levels were 
somewhat more significant (Figure 10c). Specifically, two 
wider coherent bands (9-12 days and 4.5-5 days) were observed 
with an additional peak at 3.6 days. P2 led P3 by about 579 
for the first band while P3 led P2 by 8° and 48° for the 
second band and by 1009? for the 3.6 day signal. These signals 
corresponded with the maximum energy levels observed at P2 500 


m (12 and 3.6 days) and those at P3 (9, 5 and 3.3 days). 
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Figure 10a. 
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Segment 1 Coherence and Phase for Alongshore 
Current Components: Coherence and phase are 
Shown for P4 with P2 (top) and P4 with P3 
(bottom). The solid line on the coherence plot 
indicates the 95% significance level. A 
negative phase means that the first series 
leads the second. 
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Figure 10b. 
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Segment 1 Coherence and Phase Between P2 and P3: 
Coherence and phase are shown for 350 m (top) and 
500 m (bottom). Other details are the same as in 
Figure 10b. 
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Figure 10c. Segment 1 Coherence and Phase for Across-shore 
Current Components: Details are as in Figure 10a 
except this figure shows the across-shore 
component coherence and phase between P2 500 m 
and P3 500 m only with three peaks from 9-12 
days, 4.5-5 days and 3.6 days. 

Therefore, 1t appeared that the lowest frequency across-shore 

component signal (periods between 9 and 12 days) propagated 

offshore from P2 to P3 while the higher frequencies (periods 


less than 5 days) moved onshore from P3 to P2 for both 


alongshore and across-shore components. 


B. SEGMENT 2 

Segment 2 spanned a time from 30 August 1990 through 17 
March 1991 (200 days) (Figure 2). Data were collected from 
the Pt. Sur and Monterey Bay moorings only (P2, P3, MB1, and 


MB2). MB2 was positioned about 30 km seaward of MB1 and its 
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position relative to the others is shown schematically in 
Figure 11. The Pt. Sur moorings provided data at 350 m, MB2 
at 416 m and MB1 at 500 m. The vertical coherence at mooring 
P2 during other time segments indicated that the currents at 


P2 350 m can reasonably be compared to currents at MB2 416 m. 
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sie 


Figure 11. Schematic Mooring Locations for Segment 2: 
Distances between moorings are shown as in 
Figure 4. Approximate location of the Monterey 
Canyon is also illustrated. The canyon axis is 
indicated by the darker shading. 


1. Time Domain Analysis 


In general, the mean flow during this time period was 
also directed poleward and offshore (Figure 12, Table 6). P2 
had a mean alongshore flow of 5.6 cm sec! with a maximum of 


I 


23.8 cm sec” and a minimum of -10.1 cm sec! (Table 6a). 


The across-shore component (mean -0.4 cm sec’? ) was directed 
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Figure 12. 
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Vector Velocities for Time Segment 2: Velocity 


vectors are shown in cm sec’! with poleward 


alongshore, as defined in the text, towards the 


top of the plot. Data are from P2 (top), P3 
(second), MB1 (third), and MB2 (bottom). 
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Table 6a. TIME DOMAIN STATISTICS FOR SEGMENT 2 
_(ALONGSHORE COMPONENT) 
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offshore with a maximum offshore current of -7.9 cm sec and 


a maximum onshore flow of 5.0 cm sec’? (Table 6b). For the 
first month, steady poleward flow dominated (Figure 12). 
Around 28 September, a strongly periodic flow appeared with a 
period between 20 and 30 days. This signal was not clearly 
observed at P3, located about 25 km seaward of P2, but was 
observed at MB1 and MB2. 


The flow at P3 was generally poleward and offshore 


Si 


with an alongshore component mean of 4.5 cm sec with a 


si F 


maximum of 14.5 cm sec and a minimùum of -12.0 cm sec 


(Table 6a). The across-shore component, with a mean of -2.2 


cm sec’? , was again directed offshore. The flow here had a 


=i 


strong offshore maximum of -17.5 cm sec around 25 February 


with a maximum onshore current of 7.3 cm sec ! (Table 6b). A 
weak (* 7 cm sec!) equatorward shift appeared around 01 
October and persisted for about 16 days thereafter (Figure 


2) The current shifted back to poleward around 17 October 
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Table 6b. TIME DOMAIN STATISTICS FOR SEGMENT 2 
(ACROSS -SHORE COMPONENT) 


Mooring | Meter | Mean Maximum Standard Mode 1 Mode 2 | 


(cas) (cm/s) Deviation Eigenfunction Eigenfunction 
(cm/s) (Eigenvalue = (Eigenvalue = 
22.88 (cm/s)2 for | 5.99 (cm/s)? for 
64.5 % of total 16.9% of total 
variance) variance) 





and remained so until 07 March when another equatorward event 


(* 12 cm sec 1) dominated for the remainder of the segment. 


At MB1, currents were directed equatorward and onshore 


with a mean equatorward flow of 0.2 cm sec"? 


and a larger 
onshore mean of 1.2 cm sec? (Tables 6a and 6b). Although 
weaker than the other moorings, the stick vector plot (Figure 
12) also showed an oscillatory flow with a period of about 1 
month with phase about opposite (180°) to the currents 


observed at P2. 


At MB2, alongshore flow was directed poleward with an 


1 ZEE 


alongshore mean of 2.6 cm sec” , a maximum of 17.7 cm sec 


and a minimum of -7.8 cm sec"? (Table 6a). The across-shore 


= 


flow (0.8 cm sec mean) was slightly onshore exhibiting a 


2 ll 


maximum of 13.6 cm sec and a minimum of -9.9 cm sec 
(Table 6b). The flow here appeared steadily poleward from 30 
August until 06 November. Beginning around 01 October the 


currents appeared visually coherent with the equatorward 
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bursts described at P2, but were less energetic. The stick 
vector plots provided initial indications of a coastally 
trapped, poleward propagating disturbance between P2 and the 
Monterey Bay moorings that was less obvious as a modulation of 
the mainly poleward flow at the offshore mooring, P3. 

EOF results during this time period were again 
consistent with the mean flow statistics. For the alongshore 
components, the first two modes explained 79.1% of the total 
variance. As in Segment 1, the mode 1 eigenfunction showed 
consistency at all moorings except MB1 (Figure 13, top left). 

All first mode alongshore component eigenfunctions had 
the same sign. Together with the stick vector plots, this 
information indicated a tendency toward poleward flow at all 
moorings. The second mode showed a tendency for P2 and the 
Monterey Bay moorings to oppose the flow at P3 (Figure 13, 
bottom left). This could perhaps be attributed to the low 
frequency fluctuation which was present at P2 and MB2 but not 
aE P3. 

The first two across-shore modes, explaining 80.4% of 
the total variance (Table 6a), were used to describe the flow 
structure. The mode 1 across-shore eigenfunction described 
offshore flow at the Pt. Sur moorings with P3 supplying the 
largest amount of energy. Tō the north, an opposing flow 
could be seen at MB1 and MB2 (Figure 13, top right). The 
opposing mode 1 flow between Pt. Sur and MB2 is significant 


and is discussed in greater detail in association with 
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are shown with vectors proportionally scaled 
for visual effect. 


Segment 23. Much like Segment 1, the mode 2 across-shore 
eigenfunction revealed an opposing structure between P2 and 
P3. MB2 followed the structure at P3 (Figure 13, bottom left) 
and appeared as an onshore flow also opposing P2. The flow at 
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MB1 contributed more to the total variance than it did during 
segment 1 but on the whole, it played a minor role in the 
overall current structure. 

The mode 1 alongshore amplitude time series showed the 
highest correlation with the P2 alongshore velocity component 
time series (.94). Thus P2 most nearly resembled the mode 1 
amplitude time series (Figure 14, top), showing P2 as the 
primary source for mode 1 variance. Mode 2 showed the highest 
correlation at P3 (.81) (Table 7). 


Table 7. SEGMENT 2 CORRELATION COEFFICIENTS FOR FIRST AND 
SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 
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The across-shore modal amplitude and time series 
correlation showed P3 with the highest correlation (.97) for 
mode 1 while MB2 exhibited the highest correlation for mode 2 
(.75) (Table 7). The P3 across-shore velocity components 
plotted with the mode 1 amplitudes showed this relationship 
very well (Figure 14, bottom). P2 and MB1 did not contribute 
significantly to the across-shore flow spatial structure. 

In general, the primary current structure showed a 
poleward flowing current with a mean slightly smaller than 


that for Segment 1. The second mode alongshore eigenfunction 
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Figure 14. A Comparison of the Mode 1 Velocity Amplitudes 
with the Best Correlated Individual Time Series: 
The top panel shows the P2 alongshore components 
and the bottom panel shows the across-shore 
components. 

revealed consistent flow structures at P2, MB1 and MB2 that 

opposed the P3 structure. The equatorward bursts visible in 


the stick vector plots at these moorings could explain the 


second mode structure. 
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2. Frequency Domain Analysis 
Variance conserving spectra were calculated for 
alongshore and across-shore velocity components. The 
alongshore component spectra showed energy peaks of 29.2 days 
at all moorings. The energy maximum at P2 (51.2 cm? sec ?) 
greatly exceeded those at the other moorings (between 9 and 12 


2 


cm? sec'?) 


(Figures 15a and 15b). 
The across-shore component spectra showed significant 
poscat 19.5 days for P2 and P3, 29.2 days for MB2 and 3.2 


days for MB1. These peaks were all less than those described 


2 


for the alongshore components (between 4 and 6 cm^ sec?) 
iBsigure 15b). 

The alongshore components cross-spectra, showed 
coherent signals between several mooring pairs. These 


generally occurred at periods of 58.5, 29.2 and 19.5 days 
(Figures 16a and 16b). The phase varied with each period and 
between the mooring pairs. Phase diagrams for these three 
periods suggest complex phase relationships (Figure 17). 
Arrows indicate the direction of phase propagation, e.g. P2 
leads MB2 by 37° at T = 29.2 days. No arrow indicates no 
Significant coherence. A marginally significant arrow is 
sometimes shown as a dashed line if the phase was steady and 


consistent with the other pairs. 
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Segment 2 Variance Conserving Auto-Spectra 
for P2 (top) and P3 


panels) and ACS 
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Figure 15a. 
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Figure 15b. 





Segment 
for MB 1 
panels) 





3 > 
c 
+ 
E S e 
= = = 
> 
5 
t - 
Z 2 
7 ^ 
* 
a 
-~ = - 
= D = 
= 3 = 
2 
> 
7 + 
== e T ^ m ^ ri € - ^ == 
T A ri — = 
Ce SANE C 
o ò 
e 
20 
-e 
N 
Z a 2 
— +. e 
— 2, = 
e 
2 
2 
5 
Tas c 
Z 2 
J un 
on 
= 
- GO - 
= 3 = 
a 
2 
i 
> 
m + 
-2 == 





Cos AUDE C 


2 Variance Conserving Auto-Spectra 
(top) and MB2 (bottom): ALS (left 
and ACS (right panels) are shown. 


61 


(t pa 


¿ph 


Across-shore Components 


Alongshore Components 


Figure 16a. 
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Segment 2 Coherence and Phase for Alongshore 
Current Components: Plots are shown for P2 
with MB1 (top) and P2 with MB2 (bottom). 
Other details are as in Figure 10a. 
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Figure 16b. 
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Segment 2 Coherence and Phase for Alongshore 
Current Components: Plots are shown for MB1 
with MB2 (top) and P3 with MB2. Other details 
are as in Figure 10a. 
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The 29.2 day signal (Figure 17, middle) represented 
the largest energy peak. At this period, P2 led P3 by 48% and 
MB2 by 38 consistent with a poleward propagating disturbance. 
Another interesting observation could be made by inspecting 
the phase between P2, MB2 and MB1. P2 led MB2 by around 40° 
and MB2 led MB1 by about 135%. Also, consistent with the 
stick vector plots, MB1 appeared to be 180% out of phase with 
P2. Adding the phase from P2 to MB2 with that from MB2 to 
MB1, yields a result consistent with the phase between P2 and 
MB1 and adds confidence to the results. 

The 19.5 day signal showed some characteristics that 
were consistent with the other two periods. Specifically, P2 
led both MB2 and P3 by 42° and 23°, respectively. In contrast 
to the other signals, MB2 led P3 by 24° and MB1 led P2 by 
170°. The energy at 19.5 days was less than that of 29.2 days 
but significantly higher than the 58.5 day signal. 

For the 58.5 day signal, P2 led MB2 by 42%, MB1 by 
1709, and MB2 by 23°. P3 was coherent with MB1 and led by 
136% and MB2 led MB1 by 132%, much like the 29.2 day signal. 
In general, energy seemed to propagate in a poleward direction 
but the significance of this picture is probably small since 
the energy level at 58.5 days was the lowest of all coherent 
periods. 

By knowing the phase and mooring separation distances, 
alongshore wavelengths between P2 and MB2 and across-shore 


wavelengths between P2 and P3 could be calculated. The 
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Figure 17. Segment 2 Phase Diagram: Arrows indicate the 
direction of phase propagation between moorings. 
The phase in degrees is shown for each period of 
Significant coherence. 


results for the most significant energy levels (29.2 and 19.5 
days) are summarized in Table 8. The phase error of cross- 
spectral analysis are also presented here. It is important to 


note the phase errors since this could introduce a significant 
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source of error in the wavelength calculations, which will be 
compared to theory in Chapter V. 


Table 8. OBSERVED PHASE AND WAVELENGTHS FOR LOW FREQUENCY 
OSCILLATION (SEGMENT 2): A positive phase speed 
indicates that the first cst mooring leads the second. 


Mooring Phase Phase Phase Phase Observed Separation 
Pair (o) (days) error speed wavelength Distance 
I AE (km) (km) 


T = 29.2 days (mis & MB2 aS 389.5 
(ALS) 


(ACS) 
(ACS) 


T = 19.5 days P2 & MB2 351.0 
(ALS) 


T = 19.5 days P2 & P3 122.3 
(ACS) 





Observed wavelengths were calculated by multiplying 
the separation distances by 360% and dividing by the phase 
difference. This procedure produced alongshore wavelengths of 
390 km for the 29.2 day signal and 351 km for the 19.5 day 
signal. The across-shore wavelengths were 192 km and 122 km 
for the 29.2 and 19.5 day signals, respectively at Pt. Sur. 
Between the Monterey Bay moorings, an across-shore wavelength 
of 78.8 km was calculated for the 29.2 day signal but with a 
phase speed in the opposite direction from off Pt. Sur (Table 
8). This procedure calculated the largest possible wavelength 
based on the phase differences. The phase speed was 
calculated by dividing the alongshore separation (41.1 km 


between P2 and MB2) and the across-shore separation (26 km 
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between P2 and P3, 30 km between MB1 and MB2) by the time 
lags. 

Cross spectra calculations showed little coherence 
between the across-shore velocity components. The only 
significant coherence occurred between moorings P2 and MB2 at 
a period of 19.5 days (not shown). MB2 led P2 by about 1759 
so the two moorings were generally out of phase in the across- 
shore component. A similar situation existed in Segment 3, 
where a mesoscale meandering feature dominated the currents at 


these moorings. 


C. SEGMENT 3 

Segment 3 lasted from 26 May 1991 until 18 November 1991 
(177 days) (Figure 2). Tar addition to the Pt. Sur and 
Monterey Bay moorings, this time segment included data from 
moorings D and E near the Farallon Islands. Mooring D was 
located about 110 km to the north of MB2 and mooring E was 
positioned 32 km to the north of D (Figure 18). Currents at 
P2 were obtained at depths of 350 m and 500 m but the other 
moorings only had data at 400 m. For purposes of EOF and 
spectral analyses, the currents at P2 were linearly 
interpolated to a depth of 400 m to approximately match the 
others (Figure 19). 

1. Time Domain Analysis 

Time domain statistics again showed a poleward 


alongshore component mean at all moorings. At P2, the mean 
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Figure 18. Mooring Separation Distances for Segment 3: 
Details are as in Figure 4. 


alongshore current was 6.3 cm sec’? with a maximum of 33.3 cm 
-1 


sec and a minimum of -19.0 cm sec"? (Table 9a). The 


across-shore statistics (Table 9b) showed a mean offshore flow 
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are shown as in Figure 5. Data are from P2 
(interpolated to 400 m; top), MB1 (second from 
the top), MB2 (third from the top), D (second 
from the bottom) and E (bottom). 
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ll 


ot 1.0 m sec with a maximum onshore flow of 9.8 cm sec"? 


and a maximum offshore flow of -18.6 cm sec"?. 


The statistics at MB2 were similar showing an 


E 


alongshore mean of 4.2 cm sec with a maximum of 23.2 cm 


sec”? and a minimum of -10.2 cm sec"? (Table 9a). The 


zd 


across-shore component had a mean of 0.9 cm sec with a 


21 ab 


maximum of 17.9 cm sec and a minimum of -10.9 cm sec 


(Table 9p). 


Table 9a. TIME DOMAIN STATISTICS FOR SEGMENT 3 
(ALONGSHORE COMPONENT) 


Mooring Meter | Mean Mode 1 Mode 2 
depth | (cm/s) iati Eigenfunction Eigenfunction 
(m) (Eigenvalue = (Eigenvalue = 
193.80 (cm/s)2 65.18 (cm/s)? 
for 64.7% of for 21.7% of 
total variance) total variance) 


Eran iic Rr ces cem em [ar 


PIE IA K 
> fa [io] aef ssl ee f ou M or 





* velocities interpolated between 350 m and 500 m 


From the beginning of the time segment until 27 June 
the flow was poleward at P2 400 m, with the maximum over the 
entire segment (33.3 cm sec `) occurring around 02 June. On 
27 June, the flow abruptly shifted to equatorward for four 
days and then reversed back to poleward until early August. 

From around 07 August to 07 September a relatively 
-1) 


strong (10-20 cm sec equatorward event dominated the 


currents at both P2 and MB2. The across-shore components at 
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these two moorings exhibited a remarkably strong negative 
correlation. Examination of the across-shore components of P2 
and MB2 alone (Figure 20), clearly showed onshore flow at MB2 
opposing offshore flow at P2 suggesting the presence of a 


mesoscale, rotating feature between the two moorings. 


Table 9b. TIME DOMAIN STATISTICS FOR SEGMENT 3 
(ACROSS-SHORE COM COMPONENT) 


Meter Mean Maximum Standard Mode 1 Mode 2 

depth (cm/s) (cm/s) Deviation Eigenfunction Eigenfunction 

(m) (cm/s) (Eigenvalue = (Eigenvalue = 
31.96 (cm/s)2 19.45 (cm/s)? 
for 41.3% of for 25.1% of 
total variance) total variance) 


| 
2 [ae | -no! 9.8 | me] eo 


lá [coso 
a id AI es 
rui cR 
|o two | os} ws | of «2 | oi |  coer | 
fe lo | og | wal wo] e | ow "m 


* velocities interpolated between 350 m and 500 m 





P2 currents in early August prior to the event were 
poleward but rotated clockwise slowly until 03 August. The 
current then shifted abruptly to equatorward and onshore 
briefly. Finally, it turned offshore on 07 August to begin 
the event. The current maintained an offshore, equatorward 
direction at P2 until the event ceased around 07 September 
(Figure 21). Following the event, the current shifted to 
poleward and onshore for approximately one month until an 


ESI 


equatorward burst of about -10 cm sec interrupted the flow 


on 28 September (Figure 19). A similar deviation occurred on 


Jl 


28 October. These equatorward signals appeared at about the 
same period as those observed in Segment 2 (30 days). 
Currents at MB2, directed poleward and onshore during 
the beginning of the event, did not appear to feel the 
feature's influence until 13 August. On this day the current 
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Figure 20. Across-shore Component Velocities at P2 and 
MB2: A mesoscale event commenced around 07 
August, characterized by negative correlation 
for the next 25 days. 
turned to equatorward and onshore and remained so until 07 
September (Figure 21). Flow also shifted to poleward after 
the event here and followed a similar pattern as that 
described for P2; i.e. a poleward relaxation to a near O0 cm 
sec? current on 28 September followed by an equatorward 
shift (-5 cm sec!) 30 days later. The flow at both P2 and 


MB2 during October and November appeared remarkably similar to 


m2 


that during the same months, one year earlier in Segment 2 
elgure 12). 

In general, poleward flow dominated from May until 
early August at both Farallon moorings. At mooring E though, 
the flow shifted to equatorward (max. -18 cm sec? ) on 22 
June and remained so for about 9 days thereafter. This 
occurred about five days prior to the equatorward shift 
described at P2, approximately 183 km to the south. However, 
moorings MB2 and D located between P2 and E, did not show such 
a shift so a connection between these equatorward 
interruptions does not seem likely. The shift at mooring E 
was probably due to a cold filament observed off Pt. Reyes, 
north of San Francisco Bay during this time segment (Ramp et 
al., 1992). Around 21 August, approximately two weeks after 
the start of the event at P2 and MB2, flow at moorings D and 
E weakened significantly and shifted from poleward to 
equatorward. This raised the question of whether or not a 
mesoscale event occurring some 140 km to the south could 
manifest itself into the flow pattern at D and E. With the 
limited data available, a connection between these two events 
seemed speculative but it does present an interesting 
possibility. 

At mooring D, the flow shifted back to poleward around 
27 September and remained weak (< 6 cm sec? ) for the 
remainder of the segment (Figure 19). At mooring E, a strong 


equatorward reversal occurred around 23 September, which was 
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Data are from MB2 (top) and P2 
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velocity vectors are shown 25 July 


Segment 3 Vector Velocities During Mesoscale 
through 04 September with details as in 


Event: 
Figure 5 
(bottom). 


Figure 21. 


not observed at mooring D. For the next month, the current 
at E grew steadily to a maximum of -18 cm sec. It then 
shifted back to poleward on 21 October and remained fairly 
weak and much more variable (Figure 19). 

In the empirical mode decomposition, the first two 
alongshore modes accounted for 66.4% of the variance. The 
mode 1 alongshore eigenfunction (Figure 22, top left) showed 
the flow at all moorings to have the same sign. As discussed 
previously, this could be interpreted as the poleward 
undercurrent seen in the stick vector plots. The mode 2 
alongshore eigenfunction showed a divergent flow between the 
Farallon moorings and those to the south (Figure 22, bottom 
left). This could perhaps be due to the rotating feature 
observed at P2 and MB2, which impacted the currents there but 
not farther north at moorings D and E. 

As with the other segments, the mode 1 across-shore 
eigenfunction at P2 and MB2 described an opposing flow 
structure (Figure 22, top right). Moorings D and E were not 
strongly correlated with the Pt. Sur and Monterey Bay moorings 
indicating that the mesoscale event at these southern 
locations dominated the mode 1 across-shore component 
Structure. The mesoscale feature described earlier appeared 
to influence the mode 2 across-shore current structure as 
well. The MB2 mode 2 eigenfunction had a sign opposite to P2 
(Table 9b) but the magnitude (-0.01) was much smaller than 


that at P2 (0.16) and thus does not show in Figure 22. 
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component time series 


Modal amplitudes were again correlated with alongshore 


(Table 10). 


P2 showed the strongest 


correlation with the mode 1 amplitudes having a correlation 


coefficient) cE 0.91. 


The P2 alongshore component plot and 


modal amplitude time series showed that this mooring 
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The across-shore modal 


Table 10. SEGMENT 3 CORRELATION COEFFICIENTS FOR FIRST 
AND SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 


P2 400 m 
MB1 500 m 
MB2 418 m 
D 400 m 
E 400 m 





contributed most of all the moorings to the mode 1 variance 
(Figure 23). Mooring E was most nearly correlated with the 
second mode alongshore eigenfunction (Table 10). P2 again 
showed the strongest correlation (0.92 coefficient) and the 
relationship between the two is illustrated in Figure 23. 
Much like the alongshore components, mooring E was best 
correlated with the mode 2 eigenfunction. 
2. Frequency Domain Analysis 

The variance conserving spectra (Figures 24a through 
24 c) for alongshore components showed that Moorings D, E and 
P2 peaked at 29.2 days, similar to Segment 2. MB2 experienced 
a maximum energy level at a shorter period of 14.6 days. 
Mooring E contained the most energy, 102.9 em- sec -) 22 had 


2 


less energy with a 68.4 cm? sec'^ peak. MB1 showed a 19.5 day 


2 


peak with an energy level of 11.0 cm? sec"*. D and MB2 peaked 
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Figure 23. 
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Figure 24b. 
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Mooring E 
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with energy levels of 9.0 and 17.6 cm? sec? respectively 


(Figures 24a through 24c). 

Across-shore component energy peaks occurred at 
periods ranging from 5.8 days (moorings D and MB1), 7.3 days 
(moorings E and MB2) and 29.2 days at P2 (Figures 24a through 
24C). The peak at P2 contained the maximum energy of all 
moorings (3970 cm? sec) (Figures 24a through 24c). The 
pervasiveness of the 29.2 day peak is interesting, but missing 


this time at MB2 in both the alongshore and across-shore 


components. 
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Figure 25. Segment 3 Coherence and Phase for Across-shore 
Components for P2 and MB2: The solid line on 
the coherence plot indicates the 95$ 
significance level. The two mooring locations 
are coherent but 180° out of phase. 
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Cross spectra showed the only coherent alongshore 
component signal between moorings D and E at periods from 9.7 
to 8.4 days. The energy levels at these periods were 
significantly smaller than at longer periods; thus the 
importance of this coherence is not clear. Significant 
coherence occurred between the P2 and MB2 across-shore 
components at periods from 58 to 19.5 days (Figure 25). These 
signals were =180° out of phase supporting the opposing flow 
Structure resulting from the presence of a rotating mesoscale 
feature. No other coherent signals were observed between any 


other mooring pairs. 
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V. DISCUSSION 


A. CALIFORNIA UNDERCURRENT 

Time series plots and EOF analysis consistently 
demonstrated a predominantly poleward flow during the three 
time segments at all moorings. These results were generally 
consistent with the flow of the California Undercurrent (CUC) 
as described by Chelton (1984), Wickham et. al (1987) and 
Robson (1990). Chelton (1984) computed seasonal geostrophic 
velocities and found the maximum poleward current occurred 
during December off Point Sur. Segments 1 and 2 time series, 
however, Showed current maxima during February. Chelton based 
his calculations on long term averages of all available data 
from the California Cooperative Oceanic Fisheries 
Investigations (CALCOFI). The reliability varied from month 
to month and from location to location and could possibly 
explain the difference between Chelton’s findings and those of 
Segments 1 and 2. Also, Chelton considered the entire water 
column to include currents extending to the surface while 
these results examined the 350 - 500 m depths only. The 
important point is that definite poleward maxima existed 
during winter in both Segments. 

The CUC signature during Segment 3 agreed well with 


Wickham et. al (1987) who found predominantly annual 
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varlability at Point Sur with a maximum poleward flow in May 
and June. Specifically, Segment 3 time series clearly showed 
a maximum at all locations to include Point Sur, Monterey Bay, 
and the Farallon Islands moorings during May 1991. During the 
first two segments, the flow at P3 provided the most 
consistent poleward flow and did not exhibit the periodic 


reversals to equatorward flow seen at P2. 


B. COASTAL TRAPPING 

Although the poleward flowing CUC dominated the mean 
currents, a low frequency oscillation was superimposed on the 
poleward undercurrent during Segment 2 (October 1990 to 
January 1991). This deviation was observed more clearly at P2 
than at P3 indicating the presence of a coastally trapped 
phenomenon. If a wavelike feature is assumed, the coherence 
and phase information presented earlier can be used to 
calculate phase speeds and wavelengths for the different 
energetic peaks in the spectrum (Table 8). These can be 
compared with the dispersion relations for different kinds of 
Simple coastally trapped waves to see if they provide a 
reasonable explanation for the frequencies, wavelengths and 
phase speeds observed. 

1. Rossby radius of deformation 

Significant evidence supported the importance of 

coastal-trapping. First, calculation of the internal Rossby 


Radius of deformation revealed that P2 lied close to one 
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Rossby radius of the shelf break while P3 was located outside 
of this radius. The internal Rossby radius of deformation 
represents the horizontal scale at which rotation effects 
become as important as the buoyancy effects (Gill 1982). An 
expression for this scale Can be derived following the 


procedure outlined in Gill (1982): 








a, < e (14) 
where C, = ygH, represents the wave speed of the nth mode and 
> NH ; 
H, = 1s the equivalent depth. An expression for the 


gax 


Rossby radius was determined by substituting the equivalent 


depth into the equation for c, and then substituting c, into 


(14) to get: 
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CTD casts conducted near P2 provided Brunt-Vaisala  (N?) 
profiles for October, November, December 1990 and February, 


1991. Using the standard equation: 


N= (- 


d» ka 
l% 
ele 

nje 


il 


a depth averaged value of 5.50 x 10^? rad sec l was found. 


Using the current meter mooring depth at P2 (800 m) for H and 


l 


the Coriolis parameter of 8.65 x 10^? rad sec'! produced a 
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first mode (n=1) Rossby radius of 16.2 km, and explained why 
coastal trapping influenced the flow at P2 more strongly than 


cita at P3 (Figure 26). 
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Figure 26. Schematic Cross Section Between P2 and P3: 
Diagram shows the bottom slope and the Rossby 
radius of deformation, drawn to scale at the top 
of the plot. 
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2. Critical Latitude for Coastal Trapping 
Additional support for  coastal-trapping in the 
vicinity of Point Sur came from a study by Clarke and Shi 
(1991). A simple formula produced the critical frequency at 

which waves become trapped at a given latitude: 


Q BC, 
E 2f 


O 





cos8 (16) 


where B is the latitudinal change in planetary vorticity 


h 


(2 x 10°** m* sec’*), c, = the phase speed for the n+” int 


n 


mode, f, = the Coriolis parameter (8.64 x 10"? sec"*) and 0 is 
the angle between the coastline and true north determined from 


a chart to be 35°. Using the expression: 


with N = 5.5 x 10°? rad sec’! and H = 800 m, gives a value of 


cC} = 1.40 m sec}. Substituting into (16) yields a frequency 


Of 9-33 x 108 sec Moneo ( T= 22) of 548 days. Waves 


with periods greater than this will leak offshore in the form 
of a planetary Rossby wave while smaller periods are coastally 
trapped. Both observed periods (29.2 and 19.5 days) were well 
within this maximum trapping period offering additional 


evidence to strengthen the argument that the low frequency 
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oscillations observed off Point Sur were trapped near the 
coast. 
3. Poleward Propagation 
Alongshore component phase observations from 

Segment 2 show that the oscillation propagated poleward at all 
periods (Figures 16 and 17). This supports the theory 
presented in Chapter I that coastally trapped phenomena 
propagate in a direction such that the coastal boundary is on 
the right in the northern hemisphere. In all cases except the 
19.5 day signal, the equatorward reversals propagated from P2 
towards the north (Figure 17). This evidence along with the 
Rossby radius information and critical trapping latitude 
prompted a more detailed analysis of possible coastal trapping 
mechanisms. The analysis focused on alongshore component 
phase observations only, since the majority of the variance 


was contained here. 


C. Possible Coastal Trapping Processes 
l. Internal Kelvin Waves 

As discussed in Chapter I, the presence of an IKW 
first required the assumption of a flat bottom with a vertical 
boundary at the coast. The solution to the momentum equations 
based on these assumptions produced two waves: one that 
exponentially decayed from the coast seaward and the other 
which grew exponentially from the coast seaward. The former 


represents the only physical solution and propagates poleward 
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along an eastern boundary. The dispersion relationship of 
Equation (12) facilitated a comparison between the observed 
wavelengths and theoretical wavelengths. Only the periods of 
maximum energy (29.2 and 19.5 days) were considered. Since o 
and H were well known compared with 1, the dispersion relation 
was rearranged to produce a solution for the alongshore wave 


number for each internal mode (n): 


DE 
— 0 
(205-29) 2) 
By substituting known quantities into (17), a theoretical 


value for the alongshore wave numbers and hence wavelengths 


for the two periods of interest could be obtained. In (17): 


T, - 29.2 days, o - T. = 2.49 x 10% rad sec? (T, = 19.5 


days, @, = 3.70 x 10'$ rad sec DF H = the P2 mooring depth 


(800 m) and N = 5.5 x 10^? rad sec'l, i.e., the same value 
used in the Rossby radius calculation. The results (Table 11) 
show that IKWs of these periods are expected to have much 
larger wavelengths than those observed. The wavelengths 
decrease with increasing n, but modes higher than 1 or 2 have 
not been observed in the ocean so this seems an unlikely 
explanation. Our N is, if anything an underestimate, and 
larger H makes the theoretical wavelength larger so this 


Cannot reduce the discrepancy. 
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Table 11. OBSERVED WAVELENGTHS COMPARED WITH INTERNAL 
KELVIN WAVE (IKW) AND CONTINENTAL SHELF 
WAVE (CSW) THEORY 


I Observed IKW Theoretical Wavelength CSW Theoretical 
_ Wavelength (km) (km) Wavelength (km) 


T, = 29.2 Days 





To = 19.5 Days 


Since N was well observed, we conclude that pure IKWs cannot 


explain the observed phase relationships during Segment 2. 
2. Barotropic Continental Shelf Waves 
These waves provide the opposite "end point" for 
comparison with the calculations. As in the IKW procedure, 
rearranging the CSW dispersion relationship from Chapter I 
provided a means to calculate a theoretical alongshore 


wavenumber and wavelength for a CSW: 


2T KA 


l1? = re ke +A 2 (18) 


where f is the Coriolis parameter (8.64 x 10° rad sec’ for 
36°-20' N), k is the observed across-shore wave number, A is 


1 dh 


the bottom slope parameter ae and w is the observed 


wave frequency. Using a depth of 800 m (P2 mooring depth), 
a dh of 1800 m, a dx of 52 km (distance to P3), A = 


4 


Eus x30? m'. Alongshore wavenumber calculations produced 


a wavelength of 28 km for the 29.2 day signal and 27 km for 


9H 


the 19.5 day signal (Table 11). Assuming that less than one 
complete wavelength fit between mooring locations, these waves 
are too short to explain the observations. The preliminary 
conclusion is that the waves are either 1) hybrid waves where 
both bottom slope and stratification are important or 2) waves 
which are scattered/reflected by the Monterey Canyon, which 
renders conventional dispersion relations useless. 
3. Limitations of CTW Theory 

The differences in Table 11 can be explained by 
examining the governing assumptions behind each dispersion 
relationship. First, IKW theory assumed a flat bottom with a 
vertical wall. This assumption is not valid since the 
Monterey Canyon passes directly between P2 and MB2 (Figures 1 
and 11). On the other hand, CSW theory assumed homogeneity 
with a sloping bottom. Hydrographic data in this region shows 
vertical density stratification (Robson, 1990) and weakens the 
homogeneous assumption. Therefore, it seems that the observed 
oscillations are produced by a complex interaction of both 
processes. 

Finally, while the wavelength calculations for IKW is 
non-ambiguous, it is possible in the case of shelf waves that 
more than one wavelength "fits" between moorings P2 and MB2. 
If the actual phase difference, for instance, was 360° + 38° 
= 398°, rather than 38°, then a wavelength of 37 km is 


implied, in closer agreement with the observed wavelength. 


az 


This idea is speculative, but does represent an alternate 
solution to the problem. 
4. Additional Evidence Supporting Coastal Trapping 

Vertical coherence and phase were calculated at P2 
from May 1989 through October 1990 (17 months). Current meter 
data at 100 m, 350 m and 500 m showed an upward vertical phase 
propagation consistent with results of Romea and Allen (1983) 
suggesting an IKW presence. Romea and Allen (1983) focused on 
low frequency behavior with periods similar to those observed 
during Segment 2 (T = 20 - 40 days). They employed a model 
that examined "the response of a rotating, stratified, f-plane 
ocean with a rigid lid, forced by an alongshore wind stress at 
the coast." They calculated the response for three different 
types of wind stresses: winds acting for a brief time at a 
Single point, winds acting for a longer time at a single point 
and more generally, winds acting over a larger area for a long 
time. Their model showed that bottom reflections could be 
ignored if (1) the bottom is highly dissipative, (2) if the 
waves have sufficiently low frequencies and travel a great 
distance before hitting the bottom or (3) if the waves are 
internally dissipated prior to reaching the bottom. They also 
showed that a poleward travelling wind forced a poleward 
propagating component in the form of a coastally trapped IKW 


that had a phase velocity directed upward and a group velocity 


23 


directed downward i.e., consistent with the results observed 
dt ES 

Several other models have addressed CTW forcing 
mechanisms. Chapman (1987) applied wind forced, long coastal- 
trapped wave theory in an attempt to hindcast quantities 
observed during the Coastal Ocean Dynamics Experiment (CODE) 
which took place in the vicinity of Point Arena during 1981 
and 1982. He found that bottom pressure and alongshore 
current were forced by winds between Pigeon Point, (200 km 
south of the CODE site) and the CODE site itself. Battisti 
and Hickey (1984) also found that sub surface pressure (SSP) 
response off Washington and Oregon resulted from wind forcing 
between San Francisco and Cape Mendocino over the entire low 
frequency band for periods greater than 5 days. They found 
that local wind forcing was negligible during the summers of 
1972 and 1978, and thus remote forcing was responsible for 
fluctuations in SSP observed then. Finally, Lopez and Clarke 
(1989) considered a solution forced by combined local and 
remote winds where the local solution disregarded the 
alongshore pressure gradient but allowed local winds and the 
remote solution assumed local winds were zero and the 
alongshore pressure gradient P, * O. They considered periods 
between 22 and 45 days. Although these models focused on 
continental shelf phenomena, the results consistently show 


that remote wind forcing, south of the current observation 
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resulted in coastal trapping at periods in general agreement 
with those observed during Segment 2. 

Vertical coherence and phase were also calculated at 
P3 with data from October 1990 through May 1991 (8 months). 
During this time period current meters were located at depths 
of 100 m, 350 m, 500 m and 1000 m. In contrast to the results 
at P2, a downward vertical phase propagation was observed for 
all current meter combinations. This indicates that P3 was 
more affected by surface phenomena such as meanders and 
eddies, than by coastal trapping. 

As discussed in the previous chapter, other 
Significant dynamical differences between P2 and P3 were 
observed during the first two time segments. To review, the 
current at P3 during Segment 1 showed a persistent poleward 
flow interrupted only by a strong equatorward event toward the 
end of the segment. This event appeared at the surface and 
extended down to at least 1000 m (unpublished data). 
Additionally, the equatorward reversals observed at P2 and the 
Monterey Bay moorings during Segment 2 were less apparent at 
P3, i.e., currents at P3 appeared to be influenced less by the 
coastal boundary. This makes sense in light of P3's location 
well outside of the Rossby radius. With this evidence, we 
conclude that data from P2 and P3 were representative of two 
dynamically unique environments: an upper slope or near shore 


regime (P2) and a lower slope or offshore regime (P3). 
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5. Further Possible Analysis 

A more rigorous analysis of this problem would require 
a dedicated numerical model that considers T I 
bathymetry and density stratification. Wilkin and Chapman 
(1990) used such a model to study freely-propagating CTWs. 
Their model attempted to simulate scattering of CTWs at a site 
on the East Coast of Australia. They found that CTW 
scattering strength varied with both bottom topography and 
density stratification creating an extremely complex 
interaction between these processes. The region studied in 
Segment 2 of this analysis appeared to be influenced by 
Similar complicated interactions. A model specifically 
designed to account for the highly variable bathymetry of the 
Monterey Canyon would be necessary to produce more meaningful 


results. 


D. THE MONTEREY BAY MEANDER 
1. General Description 

Satellite imagery from the Advanced Very High 
Resolution Radiometer (AVHRR) on the NOAA 11 satellite reveals 
that a distinct meander dominated the surface flow on 01 
September 1991 (Figure 27). This picture provides the best 
image during the time frame of the energetic feature which 
dominated Segment 3. The image exhibits characteristics 
Similar to those observed by Tracy (1990). In particular, a 


warm core, anti-cyclonic meander with a radius of 
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Figure 27. 
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AVHRR Sea Surface Temperature Image for 

01 September 1991: The different gray shades 
correspond to different temperatures with white 
representing the coldest water (ale C) and 
black the warmest water (875 Gime Velocity 
vectors from 400 m at P2 (to the south) and MB2 
(to the north) are shown in white and indicate 
equatorward flow at both moorings. Vectors are 
proportionally scaled for visual effect. 
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approximately 65 km can be inferred off Monterey Bay from the 
temperature gradient in the image (Figure 27). The meander is 
bounded to the south by a pronounced cold filament extending 
some 150 km ofr Point Sur and to the north by a pocket of cold 
water off Ano Nuevo probably due to intense upwelling often 
observed there. Tracy (1990) described conditions that 
allowed a similar meander to move shoreward and cause radical 
temperature changes inside the Monterey Bay. 

It is important to examine whether features such as 
this impact the currents at 400 - 500 m depth. Ramp et al. 
(1991), Huyer et al. (1991), and Kosro et al. (1991) described 
a similar meander and cold filament off Point Arena in June 
1987. The Point Arena meander displaced isopycnals to depths 
greater than 500 m. Geostrophic velocities calculated by 
Robson (1990) confirmed that surface flow reached depths of 
500 m just outside of Monterey Bay. Additional measurements 
at mooring P3 at other times show strong eddy-like features 
extending to depths > 1000 m during spring 1990 and 1991. The 
current patterns described in Segment 3 also suggested that 
the surface feature extended to the instrument depths at P2 
(S00 m) and MB2 (418 m). 

Poleward currents dominated the flow at both P2 and 
MB2 prior to 07 August (Figures 19 and 21). It is believed 
that the shoreward movement of the meander caused this 
poleward flow to reverse at P2 first and then MB2. From 07 to 


24 August, currents were onshore at MB2 and offshore at P2 


J8 


suggesting that MB2 was located in the northern boundary of 
the meander while P2 was in the southern boundary. From 14 
August until 03 September, the flow was equatorward at both 
moorings indicating that the surface jet seen in the satellite 
image extended down to the current meter depths at P2 and MB2. 
The mean current vectors from 01 September at 400 m are 
superimposed on the image, and verify this equatorward flow on 
the edge of the meander. 

Temperature records at P2 and MB2 also showed the 
meander’s influence at depth. As the meander moved onshore, 
the temperature at the two instruments increased steadily 
(Figure 28). After this rise frequent temperature 
fluctuations, perhaps associated with the highly variable 
frontal position, were observed at both moorings (Figure 28). 
As the meander retreated, both moorings experienced a sharp, 
simultaneous drop in the temperature around 24 August possibly 
due to the cold tongue of water surrounding the meander 
(Figure 27). 

The Segment 3 vector velocity plot (Figure 19) 
presented another remarkable issue. Shortly after the onset 
of the deep meander near Monterey Bay, steady poleward flow at 
the Farallon moorings D and E, weakened significantly and 
actually reversed to equatorward. Could deep interaction 
nearly 200 km to the south somehow "shut off" the source of 


the poleward flowing California Undercurrent? Although at 
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Figure 28. Temperature Time Series at P2 500 m and MB2 418 
m: Plots show the rise and subsequent drop in 
temperature at the moorings due to across-shore 
excursions of the Monterey Bay Meander. 


first glance there appeared to be some connection, it would be 


imprudent to speculate further without additional evidence. 


E. EXPLAINING THE EOF’S 
1. Mode 1 California Undercurrent 
EOF analysis provided additional evidence to support 
the presence of the CUC. Throughout all three time segments, 
first mode EOF results consistently depicted a poleward flow 
at all locations except MB1, located inshore of the CUC in the 


Monterey Canyon. The utility of this technique relied on the 
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orthogonality principle; i.e., eigenvalues of a covariance 
matrix are orthogonal to each other and allow the flow to be 
separated into distinct patterns. EOF analysis works best 
when the first few modes explain most of the variance. 
Wallace and Dickinson (1972) described this as a necessary 
criteria to establish physical significance. The first mode 
alongshore eigenvalues explained between 64% and 66% of the 
variance with higher modes decreasing significantly 
(Figure 29). Therefore, the mode 1 result is attributed to 
the dominant poleward flow of the California Undercurrent. 
2. Segment 2 

The mode 2 EOF structure is most likely due to the 
presence of significant deviations from the mean flow. In 
Segment 2, the mode 2 structure can be attributed to the 
coastal trapping phenomenon believed to dominate the segment. 
Figure 13 (bottom left) shows this idea the best. 
Specifically, the mode 2 structure at P3 opposed that at P2 
and MB2 suggesting that the trapping occurred most efficiently 
shoreward of P3. 

3. Segment 3 

In Segment 3, both mode 1 and mode 2 eigenfunctions 
displayed characteristics consistent with the Monterey Meander 
that dominated the flow at P2 and MB2 (Figure 22). In the 
alongshore component, the mode 2 eigenfunction showed that 


both MB2 and P2 had the same sign which opposed the flow at D 


101 


Figure 29. 
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Percent Variance Decrease for Alongshore 
Eigenmodes: Shown are the percentages of variance 


explained (ordinate) vs. mode number (abscissa) 
for Segment 1 (top), Segment 2 (middle) and 
Segment 3 (bottom). 
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and E. This could be attributed to the strong equatorward 
influence of the meander as described above. In the across- 
shore component, mode 1 showed this structure the best with 
flow at P2 opposing that at MB2. This can be explained by the 
onshore flow at MB2 opposing the offshore flow at P2 while the 


meander passed through these moorings. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 

Three time segments of sub-thermocline (350 - 500 m) 
current meter data along the central California coast were 
examined to determine alongshore low frequency variability on 
time scale of two days to several months. Time and frequency 
domain analyses yielded three basic conclusions: (1) the 
poleward flowing California Undercurrent (CUC) dominated the 
mean flow during all segments; (2) the presence of a coastal- 
trapped significantly influenced the flow between P2 and MB2 
during Segment 2 and (3) the deep flow at P2 and MB2 was 
influenced by a strong equatorward meander of the California 
current: 

The mean flow in all segments closely resembled that 
expected for the CUC. Segment 1 (December 89 - April 90) most 
nearly followed this mean pattern with few significant 
deviations. The first empirical mode in all segments 
consistently revealed a poleward current structure. Seasonal 
maxima occurred during February and May and showed close 
agreement with other findings in the vicinity of Point Sur. 

Segment 2 (August 90 - March 91) displayed a significant 
deviation from the mean flow in the form of a coastal-trapped 
wave. The observed wavelengths between those for an internal 
Kelvin wave and a continental shelf wave were suggestive of a 


complex interaction between density stratification and the 
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complicated bathymetry of the Monterey Canyon. The wave 
influenced P2 and MB2 much more than P3, located outside of 
the internal Rossby radius. Though not specifically addressed 
in this research, wave forcing models developed by other 
researchers yielded results not inconsistent with the possible 
forcing mechanism in the Segment 2 wave. Further study of 
coastal-trapping should include an intensive numerical 
analysis that considers both density stratification and 
complex bathymetry. 

Segment 3 (May 91 - November 91) represented another 
deviation from the mean poleward flow. An anti-cyclonic 
meander, clearly visible from satellite imagery, appeared to 
influence the flow down to 500 m. Cross-shore velocity 
component analysis at P2 and MB2 showed a strong negative 
correlation indicating the presence of an anti-cyclonic 
rotating feature. Additionally, temperature time series 
fluctuated in a manner that clearly showed the meander's 
influence. The meander appeared to disrupt the California 
Undercurrent some 200 km to the north at the Farallones 
moorings. An intensive hydrographic study encompassing 
Monterey Bay would provide insight into the surface feature- 
depth interaction and possibly improve the understanding of 


the California Undercurrent forcing. 
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